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We study a block spin transformation in the SU(3) lattice gauge theory on anisotropic lattices 
to obtain Iwasaki's renormalization group improved action for anisotropic cases. For the class of 
actions with plaquette and 1x2 rectangular terms, we determine the improvement parameters as 
functions of the anisotropy £ = a s /a t . We find that the program of improvement works well also on 
anisotropic lattices. From a study of an indicator which estimates the distance to the renormalized 
trajectory, we show that, for the range of the anisotropy £ ~ 1-4, the coupling parameters previously 
determined for isotropic lattices improve the theory considerably. 
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I. INTRODUCTION 

Improvement and anisotropy are two key ingredients of the recent developments in lattice QCD. In QCD with 
dynamical quarks, improvement of the lattice theory is essential to perform a continuum extrapolation of light hadron 
spectra within the computer power currently available [1,2]. At finite temperatures, the expected 0(4) scaling around 
the two-flavor chiral transition point is reproduced on the lattice only with improved Wilson-type quarks [3,4]. Im- 
proved actions are applied also for staggered-type quarks to reduce lattice artifacts [5,6]. However, in order to perform 
a continuum extrapolation of thermodynamic quantities, we need to further increase the temporal lattice size A^ [7,8]. 
This requires quite large spatial lattice sizes to keep the system close to the thermodynamic limit, and the task slightly 
exceeds the current limit of the computer power for QCD with dynamical quarks [9]. 

Recently, we proposed to apply anisotropic lattices for reducing the computational demand in thermal QCD [10]. 
Because the dominant part of the lattice artifacts in the equation of state (EOS) is due to the finite temporal cutoff, 
an anisotropic lattice with a larger temporal cutoff will provide us with an efficient way to calculate thermodynamic 
quantities. We tested the idea for the case of the SU(3) gauge theory with the standard one-plaquette action. From 
a series of simulations at AT t /£ = 4, 5 and 6 with the anisotropy £ = a s /a t = 2, where a s and at are the spatial and 
temporal lattice spacings, we find that the lattice artifacts in EOS are much smaller than those on the corresponding 
isotropic lattice, and the leading scaling relation is satisfied from the coarsest lattice. This enabled us to perform a well- 
controlled continuum extrapolation of EOS in QCD. Anisotropic lattices have been employed also to study transport 
coefficients and temporal correlation functions in finite temperature QCD [11-13]. In these studies, anisotropy was 
introduced to obtain more data points for temporal correlation functions. At zero temperature, anisotropic lattices 
have been employed to study charmonium states [14-16], heavy hybrids [17], glueballs [18], and also the pion scattering 
length [19]. 

Combination of the ideas of improvement and anisotropy is not straightforward, however. The main difficulty is 
the large number of coupling parameters in improved actions on anisotropic lattices. Even in the simplest case of the 
renormalization-group (RG) improved gauge action by Iwasaki [20], which contains plaquette and 1x2 rectangular 
terms only, we have 5 parameters on anisotropic lattices, instead of 2 for the isotropic case. We have to fix them as 
functions of two parameters, the gauge coupling (3 which controls the overall scale, and the anisotropy £. Because 
the redundant parameters have no physical effects in the continuum limit, they have to be determined through a 
requirement of improvement, i.e. minimizing lattice artifacts in physical observables away from the continuum limit. 

Concrete form of the dependence on the scale and anisotropy in the coupling parameters is important for a calcula- 
tion of thermodynamic quantities [21-25]. In a Symanzik-type improvement program, it is easy to see that, at the tree 
level of perturbation theory, the coupling parameters are independent of £. Accordingly, studies of finite temperature 
Symanzik-type improved actions on anisotropic lattices have been done assuming isotropic improvement parameters 
[26-28]. When we improve the theory beyond the tree level, we have to take into account the ^-dependences in 
the coupling parameters. Isotropic parameters have been adopted also in a study of the RG- improved action on 
anisotropic lattices [11], however, without justifying the choice. 

In this paper, we study the anisotropic improvement parameters for the RG- improved gauge action. Following 
Iwasaki's program of RG-improvement using a block spin transformation, we determine the values of coupling pa- 
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ramcters which minimize the lattice discretization errors. After a brief explanation of the RG-improved action in 
Sec. II, the anisotropic gauge action we study is defined in Sec. III. We then study Wilson loops under a block spin 
transformation in Sec. IV. In Sec. V, RG-improved actions on anisotropic lattices are determined and a practical 
choice of the improved action for numerical simulations is discussed. We conclude in Sec. VI. 



Various lattice actions are expected to belong to a common universality class having the same continuum limit. 
For the SU(3) gauge theory, a lattice action may contain, for example, 1x2 rectangular loops, 2x2 squares, etc., 
in addition to the conventional plaqucttes. One combination (3 — 2N c /g 2 of the coupling parameters is the relevant 
parameter which reflects the freedom of the lattice spacing, while other coupling parameters arc redundant in the 
continuum limit. The objective of improvement is to find the values of redundant parameters for which physical 
observables from a coarse lattice are closest to their continuum values. 

In order to discuss the improvement of a lattice action, we consider RG flows of a block spin transformation as 
shown in Fig. 1. The block spin transformation halves the correlation length in lattice units but does not change 
the long-range properties of the system. Then, the coupling parameter moves toward smaller (3 corresponding to the 
correlation length becoming shorter. In this figure, Ci, c 2 , etc. denote the redundant coupling parameters, and the 
points on the 1/(3 axis correspond to the standard one plaquette actions. The hyperplane at [3 — oo (g = 0), on which 
the correlation length diverges and the continuum limit can be taken, is called the critical surface. In this surface, the 
coupling parameter does not go out of the surface of (3 = oo under the block spin transformation, since the correlation 
length after the block spin transformation is also infinity, which means that an RG flow around (3 = oo runs parallel 
to the critical surface, except in the vicinity of a fixed point (FP), at which a coupling parameter does not change 
under the block spin transformation. Therefore, an RG flow can connect only at FP to the critical surface. 

Moreover, because RG flows can be regarded as lines of constant physics, the distance between each RG flow 
becomes wider as (3 increases corresponding to physical quantities becoming more insensitive to the redundant coupling 
parameters, c±, C2, etc., as the continuum limit is approached, and only one RG flow which has properties in the 
continuum limit, called the renormalized trajectory (RT), can connect to the critical surface at FP by infinite block 
spin transformations. 

The actions on the RT are "perfect actions" which reproduce continuum properties from the shortest distances on 
the lattice [29]. If an infinite number of coupling parameters are admitted, a perfect action is a goal of improvement. 
In reality, we are forced to keep the interactions as simple as possible in numerical simulations. Hence to find the 
nearest point to RT in the restricted coupling parameter space is the problem in practice. Iwasaki applied the 
program of improvement to the SXJ(N C ) gauge theory in the weak coupling limit for the case of the action with 
plaquette and 1x2 rectangular terms [20] . He found the nearest point to the fixed point in the coupling parameter 
subspace by calculating Wilson loops perturbatively on the lattice consisting of blocked link valuables after block spin 
transformations. Another approach of an RG-improved action is the classical perfect action approach in [29], which 
is suitable for increasing the coupling parameters. See [30] for a trial to include quantum corrections. A classically 
perfect action on anisotropic lattices was studied in [31]. 



On isotropic lattices, the RG-improved gauge action by Iwasaki [20], which consists of plaquettes W^ lx1 ^ and 1x2 
rectangular loops W^- lx2 \ is defined by 



where Co + 8ci = 1 for normalization and ci = —0.331 (—0.293) to optimize the action after one (two) block spin 
transformation(s) (see below for details). This action has been shown to lead to better rotation symmetry of heavy 
quark potential than the standard one plaquette action [2,32], and to suppress lattice artifacts associated with Wilson- 
type quarks at finite temperature [3] . This action was also reported to be efficient in suppressing chiral violations in 
domain-wall quarks [33]. In two-flavor full QCD with clover-improved Wilson quarks, the first systematic studies of 
the light hadron spectrum [1] and the finite temperature equation of state with Wilson-type quarks [4,8] have been 
carried out. 

The generalization of Iwasaki's action to an anisotropic lattice is given by 



II. PROGRAM FOR RG-IMPROVED ACTION 



III. GAUGE THEORY ON ANISOTROPIC LATTICE 
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where we set eg + 8cf = 1 and c + Ac\ + 4c| = 1 for normalization. This form of the tree-level Symanzik-improved 
action with mean field improvement was studied in [26-28]. Here, we study the RG transformation of this action to 
obtain an RG-improved action. 1 

Let us denote the lattice spacing in the ^-direction as a M and the lattice size as N^. We consider the case a\ = a 2 = 
as = a s and = a t , and N\ = N 2 = N$ = N s and N4 = N t with sufficiently large N s and N t . Identifying the gauge 
field by U^(x) = exp[i(/a Ai ^4 Ai (x)], the conventional gauge action is recovered in the classical continuum limit when 



Ps = 



2N C 



2N C £ 



(3) 



where £ = a s /cit- 

We perform the Fourier transformation of A^(x) by 



A„(x) = [ Q ikx+ik » a »/ 2 A^{k), 
Jk 



(4) 



where J k = -^JL_ J]J =1 £ v K = I^S and U is integer. Then, the action reads 



S = \ala t f[ E j E - C ^ + ^^ k )~m- k ) + 

M=l k^ ^i<j,a 

J> - c\a\k\ - 4afk 2 4 )f? 4 (k)f? 4 (-k) \ + 0[A% 



(5) 



with k n = (2/a^) sin(fc M a M /2) and / M „(£;) = i(kfj,A u (k) — k 1/ A /1 (k)). We adopt the lattice Lorentz gauge by adding the 
gauge fixing term: 



S gf = a 3 s a t tr 



E W*) 



A M / (x) = {/O) - C/t/( x _ a^U^x - a^)}/a^ 



(6) 
(7) 



In order to simplify the notation, we redefine lattice momenta and the gauge field absorbing the lattice spacings 
as k^o,^ > k^ 7 k^cifj, — > fc M , and A^a^ — > A^, in the following. Then, the lattice propagator, (A®(k)Al(k')) = 
S a ,bS(k + k')D^(k), is given by 



1 1 

D ^( k ) = E 7<Hi( k ) k ?h ' + tiu{k)kl5 l0 - -(q t3 (k) - l)hk 
1=1 ^ ^ 



^44 1 W=E^4(fc)fc ; 2 +e 3 fc 



(8) 



;=i 



1 Because our action contains couplings extending over two time slices, unphysical higher lying states may contaminate correla- 
tion functions at short distances comparable to the extent of the action, as observed, e.g., in a study of glueball spectrum using 
a Symanzik improved gauge action [34]. Although these unphysical states do not affect physical properties at long distances, 
a caution is required when we have to study short distance correlators to extract physical quantities. 
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where 

gy(fc) = 1 - c{(k 2 + for i ? j (= 1, 2, 3), 

qa{k) = q 4i (k) = 1 - c\k\ - 4fc| (g) 

with q^(k) = 0. 

We consider the following Wilson loops, 

x 1) = (l/N c )tr[Un(x)U„(x + MU^x + (10) 

W M „(2 x 1) = (l/iV c )tr[C^(a;)C^(a; + A)M* + 2£)I7t(a; + A + *>)t^(z + *)E/J(a;)], (11) 

W^chair) = {l/N c )tr[Un(x)U„(x + p)U p (x + A + ^(x + v + p)U}(x + &)U}(x)], (12) 

W„„„(3 dim) = {l/N c )ti[U^x)U v (x + (i)U p {x + fi + Pp^x + v + p)^(x + p)U\{x)\. (13) 

To the leading order of perturbation theory, we get [20,35] 

N 2 — 1 

(W(C))^l-g 2 ^-F(C), (14) 

4 

F^(chair) = {N^Nt)- 1 J] + D WlW (fe) 

— 2^ f y kpDvp(k) ~\~ k v k pD 'unify — kfj,k v kpD ^pite) k p ki/k p F) pu (k))], (16) 

4 

F^„(3 dim) = (A^)" 1 [] E^ 1 - ^/ 4 )^W(*0 + (! - Kl*) D »P^ k ) + (! - #/4)JW(*)]. (17) 
for SU(N C ) gauge theory, where D plJ ^ u {k) = k*D vv (k) - k p k u D vp {k) - k v k p D pv [te) + k 2 v F> - pp {te) . 

IV. BLOCK SPIN TRANSFORMATION 

The purpose of this study is to find a fixed point, at which the parameters in the action do not change. As seen 
in the previous section, F(C) are functions of the redundant parameters c s J* , hence, if a block spin transformation 
is performed in the vicinity of the fixed point, the values of F(C) should not change. In this section, we calculate 
Wilson loops on the blocked lattice after block spin transformations in the g — > limit and discuss the fixed point. 

Following Iwasaki [20], we consider a simple block spin transformation from iVes-th to (Nqs + l)-th blocking, of 
the form 

4^BS + D (n0= -L^^Bs) (n)j (18) 

nEn' 

where we block 2 4 links at the sites n — 2n' + ^2 e p jj,, (e p = 0, 1) to 1 link at n' on the blocked lattice. The lattice 
spacings change from a p to 2a p by this transformation, while the anisotropy remains the same. Note that the scale 
factor 2 is multiplied on the right hand side of Eq.(18) to scale back to the original lattice spacings, so that the 
relevant coupling g remains constant. 

Link variables on the blocked lattice are defined by U p NBS \n) = exp(iga p ^l| i JVBS ' l (n)). Wilson loops consisting of 
blocked links are given by [20] 

N 2 — 1 

(TT (jVBs) (C)) =l-g 2 c F( Nb *\C), (19) 

4iV c 
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+ (fc( ArBS )fc( JVBS )) 2 J D^ BS )(fc) - fc( 1 ArBS )(fc( AfBS )) 2 fc( JVBS )^ BS )(fc) 
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cr=l fe„ 
^(^bs)\2 



1 - 
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+ 1- 
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Cki NBs Y 



4») 
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to the leading order, where 

fc (iVBs) = ^n bs ) = 2 sin (^JV Bs) /2)j 

ATbs-1 1 4 

ff (iVBs) (fc)= TJ _Q(l + C os(2 M ^)) 



(22) 



(23) 



M=0 ju=l 

D ( ^(k) = (k^) 2 D v „(k) + (ki NBs) ) 2 D„(k) 

-2^ ArBS )fc( A ' BS )cos((2 A,BS - 1 - l/2)fc Al )cos((2 ArBS - 1 - l/2)K)D^(k). 



(24) 



The derivation ol Eqs.(20), (21) and (22) is given in Appendix A. 

In Table I, we list the numerical results of F^ Nbs ^ for the case of the standard one plaquette action. We find that 
the values of F^ 13 ^ approach to specific values in the Abs — ► oo limit. 

Wilson loops in the limit of infinite Abs can be evaluated as follows. At long distances, the gauge propagator 
should behave like 



(A-(x)AU0)} 



in physical unit. In lattice units, it reads 



4it 2 x\ + x\ + x\ + x\ 



+ 0(l/x 4 ) 



(A»^(0)) = S^SaMn) + 0(l/n 4 ), 
/i( " ) = in? + ni + n§ + e-^' ^ 



(25) 



(26) 



4tt 2 £ 2 n 2 + £ 2 n|. + t?n\ + n\ ' 



The non-leading term of the right hand side of Eq.(26) do not contribute to the expectation value in the Abs — ► oo 
limit [20]. Hence we can neglect the higher order terms. Then the resulting Wilson loops do not depend on the 
improvement parameters cj l in the original action, since the leading term does not depend on them. 
Now, the (7 x J) rectangular Wilson loops in the limit Abs 00 ar e given by 



F^(IxJ)= Jim 



1 

JVbT"oo V S Nbs 
7-1 



E 



7-1 



2If„(m -n) + 2Jf u (m - n) + 4^(7 - k)f^2 NBS kfi + m-n) 

k=i 

+4^](J - k)f„{2 NBS kO + m - n) - 27/ m (2 JVbs JO + m - n) - 2J/ i ,(2 JVBS 7/i + m-n) 
k=i 

7-1 J-l 

-4 ^(7 - k)f„(2 NBS (kfi + JO) + m - n) - 4^(J - k)f u (2 NBS (kO + Ifi) + m - n) 



k=\ 
4 fi 



fc=i 



(27) 



p=i 



= 2 Y[ / dx p (l - x p ) 



7-1 



J-l 



^(o) + + 2 E( J - fc )/p( fc A) + 2 E( J - fc )./UM 



k=l 
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.7-1 



-J/„(Ji>) - J f v {m - 2 ]T(7 - k)U(kj> + J0)-2J2(J- k)f„(kO + 7£) 



fc=l 
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(28) 
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where 



1 4 1 

/i=l,2,3(«) = J] E ( ni _ eiXl )2 + („ 2 _ e2X2 )2 + {n3 _ £32 , 3) 2 + ^2 (n4 _ ^2- ( 29 ) 

M— 1 e^ — { — 1,1} 

l 4 l 

= 4^2 II E C 2 [(ni _ £ia , l)2 + (n2 _ £2X2) 2 + („ 3 _ ^2] + (n4 _ e42 , 4 )2 • ( 30 ) 

Here, we have used limjv^oo Y^t=i /( n /2 JV ) = Jo da;/ (a;) and a relation J* dx J* dyf(x — y) = dx(l — x)[f(x) + 
f(—x)]. Similarly, we obtain, 

4 -l 

i*~> (chair) = 2 JJ / <M1 - avX/^O) + /„(0) + / P (0) - /„(£) - /„(£) - + p))> (31) 
i Jo 



cr=l 

4 



4 „1 

F^)(3 dim) = 2 JJ / dav(l - avX/^O) + /„(0) + / P (0) - + p) - / v (p + A) - /„(£ + *))• 

er=l ^° 



(32) 



We also denote in Table I. 

From the behavior that F(C)^ Nbs ' converges monotonically to F(C)(°°^ which is independent of the coupling 

sit 

parameters c/ in the original action, i.e. the starting point of RG flow, if a block spin transformation is performed 
from the point at which F(C) < - Nbs - 1 is already F(C)(°°\ we expect that the value of F(C)^ Nbs - ) does not change 
anymore, which is the property at a fixed point. Therefore we can identify the fixed point by how close the value of 
F(C) (Nbs) is to F(C)(°°). 

Notice that the property of F(C)( Nbs ) also suggests that an RG flow from every point in the critical surface flows 
into one RG flow (RT) in finite j3, on which F(C)( Nbs ) is F(CY°°\ as shown in Fig. 1, since the starting point c/ 
of the infinite block spin transformations for F(CY°°^ must be in — oo (g = 0), 

V. ANISOTROPIC RG- IMPROVED ACTION 

We search for an action which reproduces the values of Wilson loops in the JVbs — » oo limit as much as possible 
within the restricted coupling parameter space of the action, Eq.(2). For this purpose, Iwasaki considered the average 
relative deviation of Wilson loops, 



where Y^c 1S over 4 1°°P shapes up to length 6 — plaquette (10), 1x2 rectangular loop (11), chair (12), and 3- 
dimensional loop (13) — with a uniform weight w(C) = 1/4. Eq.(33) means that, when R( Nbs *> = 0.01, for example, 
the deviation of small Wilson loops from their values in the Nbs — » oo limit is about 1% after Nbs block spin 
transformations. 

On anisotropic lattices, we generalize Eq.(33) by subdividing each loop shape into orientations and adopt a uniform 
weight for each orientation. Namely, because we have 3 spatial and 3 temporal plaquette orientations, we give 
w (spatial plaquette) = w (temporal plaquette) = 1/8. For 1x2 rectangular loops, we have 6 orientations of spatial 

loops, 3 orientations of wj^ xl \ and 3 orientations of wj.\ x2 \ Therefore, we give 1/8, 1/16, and 1/16 for their weights. 

Similarly, we subdivide 12 chair and 4 3-dimcnsional loop orientations. 

Here, we should emphasize that we are trying to reproduce the values of F^°°\C) for 10 different Wilson loops by 

controlling 3 coupling parameters for £ ^ 1 (4 Wilson loops by 1 parameter for £ = 1) at the same time, which is 

a quite non-trivial trial, and the value of R( Nbs "> indicates that F( Nbs \C) does not change within the accuracy of 
r(n bs ) under block 

spin transformations, since 

j p(A r B S )(c) approaches to f(°°\C) as iV B s increases and the change of 
F( Nbs \C) is smaller than the difference. Therefore, by measuring the indicator R( Nbs \ we can check indirectly how 
"slowly" the coupling parameters flow, i.e. how the nearest point which we find in the restricted parameter space is 
close to the real fixed point in the weak coupling limit. 

In Fig. 2, we show the iVes-dependence of R( Nbs "> for £ = 1 and 2. The results from the plaquette action (open and 
filled squares) show exponential decrease with Nbs- (Results from RG-improved actions will be discussed later.) 
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We search for the minimum point of R( Nbs *> in the parameter space (cf , c* , c|) for each value of £. Figure 3 shows 
the behavior of for £ = 2 in the subspaces c* = C2 and cf = —0.31. Fig. 3(b) suggests that the region of small 
r( n bs) S p r eads in the direction of constant cj + c 2 , which we confirm also for other cases. 

To find the minimum of R( Nbs \ we solve the equations 

OjRiNus)^ _ dF (N B3){C) F (iVBs)( C) _ F (oo) {C) 

d Cl - % 2 ^ (f(~)(c))2 ~ °' (34) 

with a — {cf,c*, c\}. We iteratively solve (34) using linear approximations 

\ 2 



E 

c 



( dF^ NBa \C) \ w{C) 



\ d Cl , s ct t J (F(-)(C)) : 



:{ci - Cjo) 



E 



aF(Afes)(cr) 



c 9c « 



F(°°)(C) - F^m)^)! 



l . L 10> L 10> L 20.' 

around (cf , c* 10 , c 20 ), where dF^ s \C)/dc i can be calculated by 



(c!o ' c ' lt " c| ° ) w(C l ), (35) 



-1 



^ = -dJ^-D^ (36) 

with given by Eq. (8). We solve the equations numerically on a 128 4 lattice. We checked that the finite volume 
effects are sufficiently small for the Wilson loops in Eq.(33). 

Results for the improvement parameters which minimize R( Nbs *> are summarized in Table II and Figs. 4 and 5 for 
Nbs = 1 an d 2. We also show the results for Nbs = in Fig. 6. In these figures, cf, c\, and c\ are shown by solid, 
dashed, and dot-dashed lines. In the followings we denote the corresponding action as RG opt (jv BS ) [the RG-improved 
action with the ^-dependent optimum values of (cf,c*,c 2 ) to minimize R^ Nbs ^}. At £ = 1, we reproduce Iwasaki's 
results [20]: cf = c\ = c\ = -0.331 (-0.293) for N BS = 1 (2). 

Figure 7 shows the values of R( Nbs > from RG opt (Ar BS ) for Nbs = 0, 1 and 2 (dashed, dot dashed, dot dot dashed 
lines, respectively). We find that the values of R( Nbs *> remain small in a wide range of £, indicating that a similar 
quality of improvement is achieved by the program of RG- improvement even at £ 7^ 1. 

Here, it is worth noting that reducing the number of independent coupling parameters has a practical benefit in 
numerical simulations. In particular, because a non-trivial ^-dependence in coupling parameters makes the calculation 
of thermodynamic quantities complicated, it is attractive to adopt ^-independent improvement parameters. 

Therefore, we study R( Nbs "> at £ 7^ 1 with the improved parameters fixed to the optimum value at £ = 1, cf = c\ = 
c 2 (= —0.331 for Nbs = !)• We denote this action as RGfi xo d(Ar BS ) [an RG-improved action with cf , c*, and c\ fixed to 
the Iwasaki's value minimizing R( Nbs '> on the £ = 1 lattice]. The result of R^ for Nbs = 1 is plotted by the solid line 
in Fig. 7. We also study the case cf = c\ — c\ = ci(£) where c\ is varied to minimize R( Nbs *> at each £. The results 
for the minimum value of R^ and the corresponding optimum value of the parameter c\ are shown by dotted lines 
in Figs. 7 and 4, respectively. For both cases, R^ becomes larger as £ deviates from 1. It means that one cannot 
keep the same quality of improvement in the whole range of £ with the RGfi xe d(Ar BS ) action nor the action with the 
constraint cf — c\ — c\ = Ci(£). 

In most simulations, however, we are interested in the cases of £ w 1-4, where the values of R^ remain O(10~ 2 ). 
In the determination of an RG-improved action, the difference between Fig. 4 for Nbs = 1 an d Fig. 5 for Nbs = 2 
is a matter of taste: Both actions are equally qualified with R^ Nbs ^> < O(10~ 2 ) and the difference in the values of 
improvement parameters should be regarded as a freedom in the choice. In this respect, we find that the variations 
of improvement parameters as functions of £ are small for £ w 1-4. 

In Fig. 8, we show R^ for various actions including the standard plaquette action and the Symanzik improved 
actions. For RG-improved actions, results are shown for RG opt (i) and RGfi X0( j(i). Similar results are obtained for 
other values of Nbs, too. We find that, although a stable improvement is achieved with the RG p t (!) action for a 
wide range of £, when we restrict ourselves to the range £ w 1-4, all RG-improved actions lead to quite small values 
of R^ <O(10 -2 ), i.e. the average deviation of small Wilson loops from the Nbs = 00 limit is less than about 1% 
after one blocking. On the other hand, for the standard plaquette and Symanzik actions, typical values of R^ arc 
0.4-0.5 and 0.25-0.3, respectively. We conclude that the RGfi xe d(jv BS ) action, in which the improvement parameters 
are fixed to Iwasaki's value for £ = 1, improves the theory well at £ w 1-4. 
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VI. CONCLUSIONS 



We studied RG-improved actions for the SU(3) gauge theory on anisotropic lattices, following Iwasaki's program 
of improvement. We determined the improvement parameters as functions of the anisotropy £ for the action with 
plaquette and 1x2 rectangular terms. We found that the program of improvement works well even on anisotropic 
lattices without losing the quality of improvement if we adjust three improvement parameters (cf , c\, c|) as functions 

ofe 

Moreover, we discussed a practical choice of improved action for numerical simulations on anisotropic lattices. 
From a calculation of an indicator which estimates the distance to the renormalized trajectory near the fixed point, 
we found that keeping the improvement parameters to the values at £ = 1 leads to the distance comparable to the 
minimum distance at the optimum (cf , c\, c|), for the range of the anisotropy £ w 1-4. This means that, for the range 
£ w 1-4 where anisotropic lattices are expected to be efficient in calculating thermodynamic quantities [10], the choice 
of Iwasaki's value for improvement parameters is acceptable also for £ ^ 1, as adopted in a previous work [11]. 

As the next step, it is necessary to confirm whether good properties of the RG-improved action at £ = 1 maintain 
also at £ ^ 1 in practical simulations, but we showed that the Iwasaki's program of an RG-improved action can be 
generalized for £ ^ 1. 
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APPENDIX A: DERIVATION OF THE WILSON LOOP AFTER BLOCK SPIN TRANSFORMATION 

We derive Eqs.(20), (21) and (22) following the appendix B of Ref. [20]. Let us introduce the Fourier transformation 

Af Bs) (a .) = I ^Hkx+kj2) A (N Bs){k ^ (37) 
J k 

where the lattice spacing for A^ BS \x) is 2 Nss a. From Eqs. (4), (18) and (37), we obtain 

4"Bs>(fc) = e-^^-^^H^H^A^ik), (38) 

where 

JVbs-1 1 4 

m»n) {k)= JJ H^'Hl). (39) 
M=0 v=\ 

We define the free propagator D ( ^ Bs) for the field i( w Bs) by 

(A< N ™\k)A b J N ™\k')) = S a ^ b S(k + k')D^\k). (40) 

We obtain 

D$»*\k) = e- 4 ( 2NBS "- 1 / 2 ) fe -e l ( 2,VBS " 1 - 1 / 2 ) fe ^( JV Bs) (fc)jDAiiy(fc)j (41) 

where H( NBS \k) is given by Eq. (24). 

The expectation value of Wilson loops that can be written as 

W{C)=Y,^ v {k)D lxv {k) (42) 
for the original lattice is obtained for the A^es-th blocked lattice by 
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W (Nbs)(q = J2c^(2 NBS k)D^ BS \k) ( 43 ) 
To derive Eq. (24), we used the fact that D /xv (k) is odd in k^ and k v when \i ^ v. 
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TABLE I. Blocked Wilson loops and those N BS 


— ► oo limit. 


F s and F t are for spatial and space-time 


Wilson loops, 


respectively. N% x N t = 


128 3 x (1280 








e = i 




F(l x 1) 


F(l x 2) 




F(2 x 2) 





0.500000 


0.862251 




1.369312 


1 


0.288104 


0.517653 




0.879783 


2 


0.216234 


0.403513 




0.720860 


3 


0.194450 


0.369800 




0.674938 


4 


0.188403 


0.360256 




0.660681 


oo 


0.186476 


0.357678 




0.658761 


5 = 2 




F S {1 x 1) 


F s (l x 2) 


F s (2 x 1) 


F s (2 x 2) 





0.673095 


1.128029 


1.128029 


1.728563 


1 


0.402015 


0.701155 


0.701155 


1.145925 


2 


0.307103 


0.555452 


0.555452 


0.952321 


3 


0.277085 


0.510586 


0.510586 


0.894124 


4 


0.268545 


0.497689 


0.497689 


0.876158 


oo 


0.265709 


0.493899 


0.493899 


0.872921 


iV By 


F t (l x 2) 


F t (l x 4) 


Ft (2 x 2) 


Ft (2 x 4) 





0.556383 


0.920967 


0.995852 


1.510346 


1 


0.348216 


0.572014 


0.655826 


1.005403 


2 


0.274425 


0.454023 


0.537939 


0.837727 


3 


0.251345 


0.418526 


0.501772 


0.788069 


4 


0.244845 


0.408394 


0.491385 


0.772639 


■oo 


0.242740 


0.405580 


0.488471 


0.770227 
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TABLE II. The improvement parameters (c 1; c l5 c|) and r( Nb s) Q f the RG-improved action RG opt (jv BS ) for Nbs = 1 and 2 
at various £. 



iVBS 




c? 


A 


c| 


^(JVbs) 


1 


0.20 


-0.349 


-0.359 


-0.073 


5.40 xlO" 3 


1 


0.25 


-0.357 


-0.356 


-0.097 


5.33 xlO" 3 


1 


0.50 


-0.363 


-0.345 


-0.204 


6.10 xlO" 3 


1 


0.70 


-0.350 


-0.338 


-0.271 


7.14 xlO" 3 


1 


0.90 


-0.337 


-0.333 


-0.316 


7.67 xlO" 3 


1 


1.00 


-0.331 


-0.331 


-0.331 


7.73 xlO" 3 


1 


1.10 


-0.326 


-0.329 


-0.341 


7.68 xlO" 3 


1 


1.50 


-0.313 


-0.324 


-0.356 


7.07 xlO" 3 


1 


2.00 


-0.307 


-0.316 


-0.350 


6.52 xlO" 3 


1 


2.50 


-0.307 


-0.302 


-0.340 


6.58 xlO" 3 


1 


3.00 


-0.308 


-0.286 


-0.331 


6.96 xlO" 3 


1 


3.50 


-0.311 


-0.270 


-0.325 


7.39 xlO" 3 


1 


4.00 


-0.314 


-0.255 


-0.320 


7.77 xlO" 3 


1 


4.50 


-0.317 


-0.241 


-0.316 


8.07 xlO" 3 


1 


5.00 


-0.320 


-0.229 


-0.312 


8.30 xlO" 3 


1 


6.00 


-0.327 


-0.208 


-0.306 


8.59 xlO" 3 


1 


7.00 


-0.333 


-0.192 


-0.302 


8.73 xlO" 3 


1 


8.00 


-0.339 


-0.178 


-0.298 


8.77 xlO" 3 


2 


0.20 


-0.360 


-0.282 


-0.122 


1.34 xlO" 3 


2 


0.25 


-0.358 


-0.282 


-0.142 


1.37 xlO" 3 


2 


0.50 


-0.328 


-0.284 


-0.225 


1.87 xlO" 3 


2 


0.70 


-0.309 


-0.287 


-0.265 


2.25 xlO" 3 


2 


0.90 


-0.297 


-0.291 


-0.287 


2.42 xlO" 3 


2 


1.00 


-0.293 


-0.293 


-0.293 


2.43 xlO" 3 


2 


1.10 


-0.290 


-0.294 


-0.297 


2.42 xlO" 3 


2 


1.50 


-0.284 


-0.299 


-0.303 


2.20 xlO" 3 


2 


2.00 


-0.283 


-0.299 


-0.304 


1.91 xlO" 3 


2 


2.50 


-0.284 


-0.295 


-0.303 


1.82 xlO" 3 


2 


3.00 


-0.287 


-0.287 


-0.304 


1.91 xlO" 3 


2 


3.50 


-0.290 


-0.277 


-0.305 


2.09 xlO" 3 


2 


4.00 


-0.294 


-0.268 


-0.306 


2.30 xlO" 3 


2 


4.50 


-0.298 


-0.258 


-0.308 


2.50 xlO" 3 


2 


5.00 


-0.302 


-0.249 


-0.309 


2.68 xlO" 3 


2 


6.00 


-0.311 


-0.232 


-0.310 


2.98 xlO" 3 


2 


7.00 


-0.320 


-0.218 


-0.310 


3.20 xlO" 3 


2 


8.00 


-0.330 


-0.205 


-0.309 


3.36 xlO" 3 
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FIG. 1. RG flow and the renormalized trajectory for the SU(3) gauge theory in infinite dimensional coupling parameter 
space. 
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FIG. 2. fl^Bs) vs . 7V BS at £ = 1 and 2 for various actions. Open symbols are for £ = 1 and filled symbols are for £ = 2. 
Results from the standard one-plaquette action are shown by squares. RG opt ( n ) is the RG-improved action which minimizes 
on the lattice with the anisotropy £. RGfl xc d( n ) is an approximate RG-improved action using the values of cj* for £ = 1. 
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FIG. 4. Improvement parameters (cf,c',C2) for the RG-improved action RG opt (i) which minimizes R w at each f. The 
dotted line is the solution which minimizes when a constraint cf = c* = c| is required. 
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FIG. 7. ij^ss) as functions of £ from the RG-improved actions RG opt (jv BS ) for N&s = 0, 1 and 2. Also plotted are the results 
for determined from RGfi xec j(i) (solid line), and the minimum obtained with the constraint cf = c* = c| (dotted line). 
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FIG. 8. ii^ 1 -* for various actions as functions of £. 
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